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Considerable progress was achieved on the quantum mechanical treatment of 
electronically nonadiabatic collisions involving energy transfer and chemical reaction in 
the collision of an electronically excited atom with a molecule. 

In the first step of our work we created a new diabatic representation for the 
coupled potential energy surfaces. In previous work, 1-3 we had developed an eight- 
configuration valence bond model of the 2 A’ manifold of NaH 2 , and we used the three 
lowest surfaces for trajectory surface hopping studies of the Na(3p 2 P) quenching 
dynamics. The trajectory surface hopping calculations were in good qualitative 
agreement with results from molecular beam experiments, thus validating the potential 
surfaces and their couplings. For more accurate dynamics calculations, though, and in 
particular for quantum dynamics benchmark calculations which may be used to test 
approximate theories, it is desirable to have a different representation. On this grant, we 
developed a two-state diabatic representation designed to realistically reproduce the two 
lowest adiabatic states of the valence bond model and also to have the following three 
desirables features: (1) It is more economical to evaluate. (2) It is more portable. (3) All 
spline fits are replaced by analytic functions. The new representation consists of a set of 
two coupled diabatic potential energy surfaces plus a coupling surface. It is suitable for 
dynamics calculations on both the electronic quenching and reaction processes in 
collisions of Na(3p 2 P) with H 2 . 

The new two-state representation was obtained by a three-step process from a 
modified eight-state diatomics-in-molecules (DIM) representation of Blais and one of the 
authors. 3 The steps involved were: 

(1) diagonalization of the 8 x 8 DIM Hamiltonian and truncation at the level of 2 
adiabatic states; 

(2) numerical transformation of the two adiabatic states and their quantum 
mechanical momentum coupling, which involves derivative operators, into an 

equivalent 2x2 diabatic representation in which the coupling is a scalar 
function without momentum operators; and 

(3) analytic representation of the resulting 2x2 numerical potential function. 

Step (2) was accomplished by the transformation procedure of Baer, 4 and in step (3) we 
used an analytic function that dissociates correctly when any atom is removed infinitely 
far from the other two. Notice that step (1) involves truncation in the adiabatic 
representation, where each electronic function is a linear combination of eight valence 
bond functions. The much simpler procedure of truncating the valence bond basis would 
not yield a qualitatively correct representation. 

Accurate fitting of adiabatic potential surfaces for systems with avoided crossings 
can be very difficult because of the sudden changes in the surfaces as avoided intersections 
are traversed. These complicated shapes must be fit accurately and in a well correlated 
fashion (i.e., the sudden change in one surface near an avoided intersection must be well 
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correlated with that in the other surface and in the coupling) or spurious behavior will be 

introduced. We have instead used the adiabatic -> diabatic transformation of Baer to obtain 
diabatic surfaces and couplings which are easier to fit We have shown that the results of an 
eight-state valence bond treatment for the lowest two adiabatic states and their nonadiabatic 

vector couplings in Na-H 2 collisions may be represented efficiently in terms of a 2 x 2 
scalar potential surface matrix. 

The new 2x2 diabatic representation has significant advantages for dynamics 

calculations over the original 8x8 DIM representation and the 2 x 2 adiabatic truncation of 
that representation in that the diabatic representation is more convenient than the adiabatic 
and in that a 2 x 2 diabatic representation is more economical than an 8 x 8 representation. 

The fact that we directly fit the 2 x 2 diabatic representation rather than recomputing it by 
numerical integration every time we need a value of a potential matrix element Ujj' means it 
takes significantly less computer time to evaluate, and this is important both for trajectory 
surface hopping and also for quantum close coupling calculations of the dynamics. In 
addition, it is much easier to transform a diabatic representation into an adiabatic one than 
vice versa because the diabatic-to-adiabatic transformation is a local one whereas the 
adiabatic-to-diabatic one is nonlocal. Thus the fact that the new representation is originally 
defined in the diabatic form allows it to be used conveniently in dynamical formulations 
based on either diabatic or adiabatic representations. We expect that it will provide a useful 
system for benchmarking dynamical methods since a variety of dynamical methods, both 
those based on diabatic representations and those based on adiabatic representations, can be 
applied conveniently with the same set of interactions. 

In addition to the specific interest in this particular system, the technique should be 
of general usefulness for representing other electronic structure calculations based on 
multiconfiguration wave functions where the number of states can be truncated to a 
reasonably small number (e.g., 2 or 3) only in the adiabatic representation, but the adiabatic 
surfaces and nonadiabatic couplings are hard to tit. 

The work just described was published in the following reference: 

“A New Diabatic Representation of the Coupled Potential Energy Surfaces for 

Na(3p 2p) + H 2 -» Na(3s 2S) + H 2 or NaH + H,” P. Halvick and D. G. 

Truhlar, Journal of Chemical Physics 96, 2895-2909 (1992). 

The second step of this project required the development of new dynamical 
methods. We developed a formalism for treating reactions with very general basis 
functions including electronically excited states. Our formalism is based on the 
generalized Newton,^ scattered wave 6-^ and outgoing wave 7 *^ variational principles that 
we have used previously for reactive collisions on a single potential energy surface, and it 
incorporates three new features. 

(i) The basis functions include electronic degrees of freedom, as required to 
treat reactions involving electronic excitation and two or more coupled 
potential energy surfaces. 
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(ii) The primitive electronic basis is assumed to be diabatic, and we do not 
assume that it diagonailzes the electronic Hamiltonian even asymptotically. 

(iii) Contracted basis functions for vibrational-rotational-orbital degrees of 
freedom are included in a very general way, similar to previous 
prescriptions for locally adiabatic functions 1 ® in various quantum scattering 
algorithms. 

We also allow for contracted translational basis functions as treated previously. 11-13 

This formalism is presented in more detail in: 

“Linear Algebraic Formulation of Reactive Scattering with General Basis 
Functions,” G. J. Tawa, S. L. Mielke, D. G. Truhlar, and D. W. Schwenke, in 
Advances in Molecular Vibrations and Collision Dynamics: Quantum Reactive 
Scattering, Vol. 2B, edited by J. M. Bowman (JAI Press, Greenwich, 
Connecticut, 1993), in press. 
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